###############################################################################################
# Causal mediation analysis: The civic effect of CRSV is not mediated by posttraumatic
# growth (PTG).
# 
# Figure I.4
# 
###############################################################################################

# DRC

par(mar=c(3,4,1,3), mfrow=c(2, 1))
plot( c(mean(((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), mean(((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2)), c(1,2), pch="", axes=F, ylim=c(0.5, 2.5), xlim=c(-1, 1), ylab="", xlab="Effect on Political Participation Index", main="DRC")
axis(1, col="white")
#grid(col="lightgrey", lwd=.5, lty=1)
abline(v=0, lty=2, col="grey", lwd=1.5)
axis(2, at=c(1:2), labels=c("ACME", "ANDE"), col="white")

denstrip((((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2), horiz=TRUE, at=2, width=.03)
denstrip((((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), horiz=TRUE, at=1, width=.03)
points(c(mean(((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), mean(((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2)), c(1,2), pch=19)

text(mean(((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2), 2, round(mean(((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2), 2), pos=1, cex=.8)
text(mean(((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2), 2-.1, paste("[", round(quantile((((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2), c(.025)), 2), " ", round(quantile((((Y.11.drc-Y.10.drc)+(Y.01.drc-Y.00.drc))/2), c(.975)),2), "]",sep="") , pos=1, cex=.8)

text(mean(((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), 1, round(mean(((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), 2), pos=1, cex=.8)
text(mean(((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), 1-.1, paste("[", round(quantile((((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), c(.025)), 2), " ", round(quantile((((Y.11.drc-Y.01.drc)+(Y.10.drc-Y.00.drc))/2), c(.975)),2), "]",sep="") , pos=1, cex=.8)
box(col="darkgrey")

# Sri Lanka

par(mar=c(3,4,1,3))
plot( c(mean(((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), mean(((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2)), c(1,2), pch="", axes=F, ylim=c(0.5, 2.5), xlim=c(-1, 1), ylab="", xlab="Effect on Political Participation Index", main="Sri Lanka")
axis(1, col="white")
#grid(col="lightgrey", lwd=.5, lty=1)
abline(v=0, lty=2, col="grey", lwd=1.5)
axis(2, at=c(1:2), labels=c("ACME", "ANDE"), col="white")

denstrip((((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2), horiz=TRUE, at=2, width=.03)
denstrip((((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), horiz=TRUE, at=1, width=.03)
points(c(mean(((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), mean(((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2)), c(1,2), pch=19)

text(mean(((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2), 2, round(mean(((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2), 2), pos=1, cex=.8)
text(mean(((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2), 2-.1, paste("[", round(quantile((((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2), c(.025)), 2), " ", round(quantile((((Y.11.sri-Y.10.sri)+(Y.01.sri-Y.00.sri))/2), c(.975)),2), "]",sep="") , pos=1, cex=.8)

text(mean(((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), 1, round(mean(((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), 2), pos=1, cex=.8)
text(mean(((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), 1-.1, paste("[", round(quantile((((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), c(.025)), 2), " ", round(quantile((((Y.11.sri-Y.01.sri)+(Y.10.sri-Y.00.sri))/2), c(.975)),2), "]",sep="") , pos=1, cex=.8)
box(col="darkgrey")
################################################################################################
